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Abstract 

A strongly correlated electron system with controlled hopping, in the line of the re- 
cently proposed generalized Hubbard models as candidates for high T c -superconductors, 
is considered. The model along with a whole class of such systems are shown to be 
completely integrable with explicit quantum i?-matrices and the Lax operators. Inspite 
of novelties in the Bethe ansatz solution, the final results do not deviate much from those 
of the standard Hubbard model. However, the symmetry of the model is changed to a 
recently discovered twisted Yangian symmetry. 
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Some unusual but universal behavior of interacting fermion systems, especially in two 
dimensions, named as the Luttinger liquid theory was speculated to be the basis for the high 
T c -superconductivity |], ||. Most significant among such properties is the separation of charge 
and spin degrees of freedom. Thus the charges of electrons are given to the pseudoparticle 
modes holons and anti-holons with charges but without spins, while the spins are given to 
the spinons with spins but without any charge. These are many-body collective modes with 
strong nonperturbative nature and unlike the quasiparticles of usual Fermi liquid which in the 
limit of vanishing interactions map into free electrons, they can not exist without many-body 
interactions. Such characteristic behavior of Luttinger liquid however is more common in one 
dimension and can be observed explicitly in Bethe ansatz solvable correlated electron models 
like Hubbard model [^, |4| . Due to this fact the investigation of correlated electron systems in 
one-dimension has become immensely important |3], |], |5], || . 

Though the study of the standard Hubbard model itself was identified to be promising |7j , 
more general models including higher nonlinear interactions with correlated hopping are being 
proposed |^, [TD[ for better description of the cuprate superconductors. Though ignoring 



certain terms and restricting coupling parameters, some exact results were obtained in || and 
[10|| , such generalized Hubbard models are not Bethe ansatz solvable and clearly do not show 
complete integrability in one dimension. 

Our aim is to propose a related though different strongly correlated electron model given 
by the Hamiltonian 

H v = - C j{ _ a) C j+1 ^ a ){t AA + if A B ~ t AA){nj{cy) + n j+l{<r)) + (tAA + t BB ~ 2t AB) n j(o-)n j+1 ( a )} 

+ Uj2 n i(+) n i(-) +h.C, (1) 



where cj, )(cj( CT )) is the creation (annihilation) operator for an electron with spin a = ± at 
site j and = cL sC^) is the number operator at site j with N a being the total number of 
sites. The coupling constants are taken as t BB = (t^ B ) 2 = e ±l2ri and scaled by putting t AA = 1, 
where the real parameters r\ and U are kept arbitrary. Notice that apart from the standard 
Hubbard interaction represented by the double occupancy term Hjj = UJ2j n j(+) n j(-) an d the 
free hopping, the model contains additional interacting terms influencing hopping. These 
are the Hirsch-like interaction || given by the second term and higher nonlinear interaction 
involving different sites given by the third term. In particular, the hopping of up (down) spin 
electron is controlled by the presence of down (up) spin electron at the same or at neighboring 
sites. Moreover, since the coupling constants are different for up and down spins as well as 
for left and right hoppings, their hopping rates might be different. It may be mentioned that 
the model proposed in JED) is of the form ([j]) with coupling constants independent of spins: 

tAB = ^ABi t%B = Ibb- 

We show that, as opposed to the model of [|1(J the generalized Hubbard model (fj) is 
exactly solvable by the Bethe ansatz with unusual features. Moreover, the corresponding Bethe 
ansatz results do not deviate much from those of the standard Hubbard model and proceeding 
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analogously one can expect to find similar Luttinger liquid like behavior. However, the present 
model exhibits a new type of symmetry discovered very recently, namely the twisted Yangian 
symmetry at the infinite chain limit. Furthermore, along with a whole class it belong to the 
completely integrable quantum systems with sufficient number of independent conservation 
laws. The associated i?-matrix and the Lax operators satisfying the quantum Yang-Baxter 
equation can be extracted in the explicit form, which are found to be intimately related with 
those of the Hubbard model. 

Using the fermionic property (n-^-b)) 2 = %(±), it is convenient to rewrite Hamiltonian ([!]) 
with our choice of parameters in the form 

H v = ~ £ 4+i(+) c i(+) ei ' ?K ' ( - )+ ^ +1C " )1 + S t +1( _ )Cj( -)e- ir?[nj(+)+nj+1(+)1 + Un j(+) n jH + h.c. (2) 
j 

We show below that the eigenvalue problem of the model can be solved exactly by using the 



coordinate formulation of the Bethe ansatz. Though this method proposed first by Bethe [11 



has been applied since to a number of models and by now has become almost an algorithmized 
problem, the present model shows surprises in its Bethe ansatz solution. 

This subtle feature is manifested already in the simple two-particle case, if we consider the 
wave function involving one up and one down-spin electrons created respectively at the sites X\ 
and x 2 . We have to distinguish naturally between the state ift( + '(xx, x 2 ) in the sector x\ < x 2 , 
i.e when the up-spin is created left to the down-spin and the state ip(~ + '(xi, x 2 ) in x\ > x 2 with 
the up-spin to the right of the down-spin. Since the interactions are short-ranged appearing 
only for the opposite spins occupying the same or the nearest-neighbor sites, all interactions 
naturally vanish when spins are placed well apart. 

Therefore for Xi,x 2 far apart the wave functions ip( ±T \xi, x 2 ) should satisfy the free 
discrete Schrodinger equation 

^ {±T \ Xl - l,x 2 ) + ij ( - ±T) (x 1 + l,x 2 ) + 4) {±T \x u x 2 - 1) + i) i±T) (x 1 ,x 2 + 1) 

= -E^\x u x 2 ). (3) 

On the other hand for the nearest neighbor occupation of the opposite spins the interactions of 
the hopping terms come into play, resulting for x 2 = Xi + 1, Xi = x, the Schrodinger equation 

^+-\x-l,x + l)+^ + -\x + l,x + l)e- iv + ifj {+ - ) (x,x)e- iv + ifj {+ -\x,x + 2) 

= -E^ + ~\x,x + l). (4) 

and similarly for x 2 = X\ — 1, Xi = x, related to ip(~ + \ Now as the Bethe ansatz demands, the 
same solution ip^ + ~~> must hold for both these equations, i.e. eqn. ([!]) at the limit x 2 — > X\ + 1 
should formally coincide with ((4|). This leads to the unexpected result 

^ (+-) (xi,x 2 ) \ XX ^X2= e~ lT1 ii(x2,x 2 ), and similarly ^ ( ~ +) (xi, x 2 ) \ Xl ^x 2 = e iv il)(x 2 , x 2 ), (5) 

where ip(x, x) = ip( + ~\x, x) = ^~ + \x, x) at the coinciding points. This unusual anyonic type 
feature introduced through coupling constant r\ can be understood also at the operator level, 
as will be demonstrated below. 
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We notice further that at X2 = xi = x both interactions involving hopping terms with 
parameter rj as well as the Coulomb term with coefficient U become active resulting the 
Schrodinger equation 

ifj {+ ~\x - 1, x)e ir > + ifj { - + \x + 1, x)e~ iri + ^~ + \x, x - l)e~ iri + ^ {+ ~\x, x + l)e ir] 

+Uijj(x,x) = —Eip(x,x). (6) 

Demanding again that (|f) should be compatible with @ when x\ — > X2 and using the relations 
(U) at the coinciding points we arrive at the consistency condition 

^ + ~\x + 1, x) - ip { - + \x + 1, x)e~ 2iri + ^ {+ -\x, x-1)- i) { - + \x, x - l)e~ 2iv 

- U^ + '\x u x 2 ) \ Xl -+ X2=x = 0, (7) 

Defining now the wave functions in the standard Bethe ansatz form 

^\x 1} x 2 ) = A^/^+p^ - A^e'^+P^, (8) 

which is valid basically in the sectors X\ ^ x 2 , one can easily calculate the energy eigenvalue 
from (^|) as E = — 2(cospi + cosp 2 )- To evaluate (H) at the coinciding points, one has to 
consider (|) yielding ^^'{x^xz) \ Xl ^ m — e~ 2lv ^~ + '(xi,X2) \ Xl -^ X2 , which in turn leads to 
the nonstandard relation 

4+~ — A + ~ = p~ 2ir >( A~ + — A~ + \ (Q\ 

P1P2 ^>2Pl C V^PlP2 ^VlPxl- V) 

Inserting ansatz (Bp in the consistency condition (|7|) and using relation (Kp one gets finally the 
two-particle scattering matrix 

S(\« A °)- (A °~ A ° )S + ^ (10) 
S{Xl ~ Aa) " A?-A^ + zf (10) 

defined as A a p \ pi = £ Cfi S(A° - X° 2 ) b c a d Af ip2 , where the obvious relations like = A±± 2 

are imposed. Here A° = sinj>j, P = | (/ + a <8> <?) is the permutation operator and 

S = diag{l,e 2i \e- 2i \l). (11) 

It is crucial to note that the S'-matrix (|10|) satisfies the well known Yang-Baxter equation 



m 

>5l2(Ai — A2)<S'i3(Ai — \3)S 2 3(\ 2 — A3) = 623^2 — A3)5i3(Ai — A3)<Si2(Ai — A2), (12) 

representing the factorizability condition for the many-particle scattering into the two-particle 
ones. This enables us to solve the general iV-body problem through the Bethe ansatz, when the 
state corresponds to the presence of total N number of electrons with M down-spin electrons. 
Each particle scattering through the rest and returning to its original position would generate 
a string of S'-matrices (|l(J) in the factorized form. On the other hand, for a closed chain with 
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the periodic boundary condition this would be equivalent to the shift operator over the total 
number of lattice sites N a . Diagonalizing this relation, as seen from fllCf ) one gets 



e 



iPjN a _ -i2r,M TT 1 "2 



n J v j 4r JL > (13) 



along with the relations 



n \ _ \0 i aU m \ __ \ j_-U 
-i2 V N 77 Aa A j 1 2 _ 77 A » A P "t" 1 2 



n A \ 2 = n r 7 ( i4 ) 

j = l A Q — ( g =1 A Q — A^ — I 



Taking logarithm of these Bethe equations we get 
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PjN a = 27rJ i + 2^tan- 1 ( — (A a -sinp j ))-2r]M 

a=l U 

N 4 M 2 

2^tan- 1 (-(A a -sinp i )) = 2ttJ q + 2 £ tan- 1 (-(A Q - Xp))+2 V N (15) 

j=l ^ /3=1 U 

where Ij , J Q are integers or half odd integers. 

We define charge and spin rapidities as Pj(qj),j = 1,2, ... ,N and X a (p a ), a = 1,2, ... , M, 
respectively, where qj = j^-Ij and p a = jf-J a are given through two independent sets of 
quantum numbers Ij and J a . Comparing (p~5| ) with the Bethe ansatz results of the Hubbard 
model we may conclude that, though due to the inclusion of various interacting terms 



influencing hopping our generalization differs considerably from the standard Hubbard model, 
the final Bethe ansatz results do not show significant changes. Nevertheless, though the 
energy eigenvalue is obtained in the same form, the determining equations for p(q) and A(p) 
are modified due to the appearance of additional phases involving coupling constant rj and 
particle numbers N and M. Going to the thermodynamic limit, analogous to the Hubbard 



model H f|, [13]], one expects to show that all low lying excitation modes of the system are 
expressed through the decoupled charge and spin degrees of freedom. 

Interestingly, for our choice of the coupling constants, the model ([!]) not only becomes 
Bethe ansatz solvable, as we have seen above, but also turns out to be a completely integrable 
quantum system with higher conservation laws like the original Hubbard model |T5|]. The 
i? Hub (A, /Lt)-matrix and the Lax operator L Hub (A) of the Hubbard model as a solution of the 
quantum Yang-Baxter equation [16[ 



Ru(Xi, X 2 )L lj (X 1 )L 2j (X 2 ) = L 2j (X 2 )L lj (X 1 )R 12 (X 1 , A 2 ) (16) 
were given in a convenient form in |L7|] as 

^ b (Aa) = {Lj\X a ) ® L^(X a ))exp(h a af +)a al )a ) (17) 

and 

^? 2 ub (A 1 ,A 2 ) = [cos 9 12 cosh h l2 {L^ + \e l2 )®L^-\6 l2 )) 

\+)l u (-)2J 



+ cos^ 12 sinh/ il2 (Li< +) (^ 12 ) ® Llt ) 0i2)){^f+)i^( Lyj)] 
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where LF {± ^ correspond to the 6-vertex free-fermionic model, h\2 = hi — /i 2 and 0u = Ai — 
A 2 , 0i2 = Ai + A 2 represent the dependence on the difference and the sum of the spectral 
parameters. The Pauli matrices <7(±) correspond to the spin up/down fermion operators q±) 
and the parameters h a are defined as sinh2/i a = jsin2A a [|17]]. The Lax-operator and the 
f?-matrix for the present model are linked intimately with those of the Hubbard model and 
can be obtained easily from flT7| ) and flTSP through twisting transformation: 

Rniv, Ai, A 2 ) = F 12 (7 7 )i?f 2 ub (A 1 , A 2 )F 12 (t?); L aj (rj, X 1 ) = F aj (rj)L^ h (X)F aj ( V ), (19) 

where the twist operator is given by 

F aj ( v ) = e^t-^wr^Ujl (20) 



Recall that the twisting transformation [jM] generates new i?-matrix and L-operator solutions 
exploiting a nontrivial symmetry of the Yang-Baxter equation. 

Another important issue is to study the effect of the additional terms is changing the sym- 
metry of this model. It is well known that the Hubbard model exhibits a Yangian symmetry 
in the infinite chain limit |1| [L9|]. It is therefore intriguing to ask whether the present gener- 
alization destroys the original symmetry completely or deforms it to another one. The origin 
of the Yangian symmetry in the Hubbard model is the rational f?-matrix of the XXX spin 
chain embedded in it, which gets associated with the algebra of the monodromy matrix at 
the infinite interval |[L9|| . The corresponding rational i?-matrix here is its twisted version (|10l), 
which gives the present model a new type of twisted Yangian symmetry, discovered recently 

ja/3 

PH , |2i]| . The expansion of the monodromy matrix: T alS = T aa 5 a /3 + Y^= aSt 1,2 
yields generators r aa , t^f of the infinite dimensional twisted Yangian algebra Y^((// 2 ), defining 



relations of which are given in explicit form in [20]. We can find a realization of this algebra 



through the fermionic operators by expressing first the generators as 

j 

fn) = T aa tT(n) (21) 



with the notations rj*£ = Uk>j(k<j) T k and T ° a = Uk T k a - Note that for T k° = 1 the 

,t 



expressions ( |2TD reduce to the undeformed Y(s/ 2 ) Yangian generators t?R = J2j tfL) - Therefore 



using the well known Yangian representation for the Hubbard model involving cL ± \, Cj(±) 
and %(±) = Cj(±-\Cj(iz), along with additional expressions 

_ 2iJ7n j( _) 22 _ -2ir)n jM j f ll _ y.22 _ /oo\ 

we can obtain an exact representation of the twisted Yangian in fermion operators. In the 



line of ||18| ) it can be shown now by direct check that the Hamiltonian (|2]) for infinite chain 
commutes with the generators of this twisted algebra along with a complementary set of such 
generators obtained by replacing 

c jH - c i( _ )e - 2 "» : , c j{+) ^{-iyc) {+v with V - -U, rj - -77, (23) 
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since under (|23D the Hamiltonian remains invariant. This proves a novel twisted Y v (gl 2 ) © 
Y_ v (gl 2 ) Yangian symmetry for our generalized Hubbard model H v (0) at the infinite chain 
limit. It is easy to check that at i] = one recovers the result related to the standard Hubbard 
model |g . 

Finally we mention about a possibility of extending the present model to an one-parameter 
family of integrable models by introducing an additional coupling constant g in the form 

11 W Z^V j+l(+) T °j+l(-) L 'j(-) c > 

3 

+ E/>ij( + )7ij(_) + h.c, (24) 
which reduces to f/,, (0) at g = and generates at g = ±77 new type of models like 

= 4+i(+) c ^+) ei2 ' 7nj ' +1C " ) + 4+i(-) c ^(-) e " i2 ^ C+) + + h -c- (25) 



Note that these models are different from those proposed in or [p7 |. Remarkably, though 
the coupling constants in ( pip are given as combinations of rj±g, the eigenvalues as well as the 
scattering matrix remain independent of the parameter g. This family of models producing 
the same Bethe ansatz results and sharing the same symmetry can be represented by the 
i?-matrices and the Lax operators, which may be obtained from those with g = through a 
simple gauge transformation 

Ri2(v,9, Ai, A 2 ) = A 12 (g)R 12 (r], X h X 2 )A^ 2 (g); L aj (r), g, Ai) = A aj (g)L aj (r], X)A a j(g) (26) 

as an operator dependent similarity transformation through A a j(g) = e^ CT (-> aCT (+) a+fT (+)j <T( -)j' ) . 

The present study of the generalized Hubbard model shows that additional terms influenc- 
ing hopping with resemblance to the recently proposed models for high T c superconductors, 
may be introduced retaining its integrability in one dimension with nearest-neighbor interac- 
tions. Moreover, for a change in operators as 

4(±) "> = ^ Wfi !w. c i(±) ^ £ i(±) = ^ w ^(±), (27) 

the generalized model (0) can be transformed back to the original form of the Hubbard model. 
However the resultant operators would show not free fermionic but anyonic type of commu- 
tation relations 

4(+)SV) + e2 %-A( + ) = °> ^ = e±2 ^ W (28) 

etc. We have witnessed the reflection of this intriguing feature in the Bethe ansatz procedure 
with the wave functions suffering phase jumps at the boundaries of two different sectors. 
Usually in the Bethe ansatz, the matching of the wave functions at the boundaries is assumed. 
However, the present model suggests for a more careful comparison with the consideration of 
their phase factors. 
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Inspite of many nontrivialities involved in the generalized model, the final Bethe ansatz 
results are not much different from the original Hubbard model, apart from a symmetry change 
and modification in the Bethe equations. This effect is like putting the corresponding vertex 
model in vertical and horizontal electric fields, which spoils the invariance of the Boltzmann 
weights under inversion of all arrows, as happens in the asymmetric 6- vertex model [22] or the 
Dzyaloshinsky-Moriya interaction 



Only in the present model the fields are not external 
but caused by the interaction of different kinds of bonds. The original Yangian symmetry 
of the Hubbard model is deformed to the twisted Yangian symmetry, while the quantum 
i?-matrix, Lax operators etc. are related through twisting transformation. 

Nevertheless, like the twisted Heisenberg spin chain [25[] or the Hubbard model with 
Ahronov-Bhom period |27| such interactions, as shown in a recent work J26[, can also be 



such interactions, as shown in a recent work 
absorbed in the boundary conditions, though in a more involved way. This also makes the in- 
vestigation of some important problems, like the influence on the effective period [P7| , changes 
in the finite temperature behavior p3 and related conformal properties 



2]§ |3 and also 

the modification of correlation functions worth studying. Similar idea can also be used for 
generating integrable coupled anisotropic spin chains |31|. Following a different approach an 
integrable coupled spin chain and a quasi two dimensional extension of the Hubbard model 



had been obtained in some earlier works [32, 33 



The author thanks Dr. Indrani Bose for valuable discussions and the referee for his con- 
structive comments. 
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